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TECHNICAL NOTE 25*^0 

APPLICATION OF RESPONSE FUNCTION TO CAICULATION 
OF FLUTTER CHARACTERISTICS OF A WING 
carrying CONCENTRATED MASSES 
By H- Serbia aod E. L. Costllow 


SUMMARY 


Concepts Involved in the baxmonlc-response-function method, such 
as the direct or conjugate characteristic modes, are illustrated by 
application of the method to the calculation of the change in flutter 
characteristics of a •win g due to adding concentrated masses. The main 
pvirpose of the numerical procedures ■which are gi'ven is to illustra'te 
the scope of the method and to make some of its abstract formulations 
more specific, rather than to stress the im mediate and current practical 
usefulness. The rigid body is' first idealized as a point mass, then as 
a distributed mass. An appendix is given which contains some of the 
essential theoretical backgroiind. 


INTRODUCTION 


In recent years, a number of papers dealing with certain special 
phases of aeroelasticlty, such as flutter, divergence, and control 
reversal problems, have been ptiblished. Each such problem has usually 
been treated as a speciaO. case so that there exist today a number of 
procedures which have little relation one to another. 

There has been in existence, however, a general approach (refer- 
ence l) to aeroelastic phenomena. Althou^ the underlying theory is 
couched in abstract terms, the method has the ad'vantages of conceptual 
sin^jllcity and physical interpretation; the abstract formulation, in 
fact, permits application to a number of apparently different problems. 

Engineers who have ■worked with ■vibration mountings are familiar 
■with the significance of the natural frequencies and modes of an elastic 
system in reference to forced oscillations. The response of the system 
•fco an external force cannot be properly appreciated ■without relating 
characteristics of that force ,to the stability characteristics (reso- 
■nances) of the system. The method described below is simply the logical 
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extension of tMs framework of ideas to the case of nonconservative 
forces as encountered in aerodynamic phenomena. 

The method of the liarmonic response function is here applied to a 
special problem, namely, the effect on flutter of the attachment of a 
rigid mass to the airplane . Specifically, it is supposed that a flutter 
analysis has been carried out on the original airplane in n degrees 
of freedom. A mass is then attached. The problem is: What are the 

flutter characteristics of the modified airplane? 

The problem is ingjortant in airplane design because (a) mass bal- 
ances are frequently attached to a control surface for flutter preven- 
tion and (b) power -pi ant changes result in the addition of concentrated 
weights. 

The results of reference 2 showed that in some instances a Hayleigh 
type analysis may show no flutter when a mass is attached, whereas a more 
rigorous analysis (based on the continuous system) would show flutter. 

The problem raised, as to how many degrees of freedom are required for 
a Rayleigh type of analysis, is one of the difficult, unsolved problems 
that do not fall within the scope of this report. 

A solution of the problem by means of the response function depends 
on the accuracy with which the system is described by the n degrees^ 
in fact, the method of the present paper is exact to the same extent as 
the n degrees of freedom describe the deflections of the system 
with and without mass. 

In the case where there are an infinite number of degrees of free- 
dom (as is the case for a continuous structure), the method of Rayleigh 
approximates the motion of the system by restricting the motion to a 
finite number of degrees. If these finite degrees are chosen to favor 
the representation of the system without mass, then the representation 
of the system with mass will be less favorable because the mass will 
Introduce discontinuities in shear which normally would not appear in 
any of the approximating modes. 

Therefore, for continuous systems one may expect that a represen- 
tation of the original system by a finite number of degrees will not 
serve so well for the system with a concentrated mass. However, the 
present analysis shows that the results of this approximation yield as 
favorable results as have been obtained by the more conventional methods 
of analysis . 

This work was done at Purdue University under the sponsorship and 
with the financial assistance of the National Advisory Committee for 
Aeronautics. 
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SYMBOLS 

abscissa of •wing elastic axis measured from wing midcbord 
point, positive rearward, in units of b 

real part of Xj 

coefficients of inertia-plus -aerodynamic virtual, work 
intensity factor ( see text) 

y~ yrS-s^s 

ris 

semichord of wing 
imaginary part of 

coefficients of elastic vlrttial work 

y ^r^rs^s 

r,s 

intensity factor (see tex±) 

N 

constant 

distance between center of gravity of attached mass and 
elastic ajd-S, positive rearward 

direct characteristic mode (normalized) 

conj-ugate characteristic mode (normalized) 

bending deflection of wing 

y 

torsional deflection of wing 
/ 

amplitude of oscUJLating force 
direct free mode 
conjiogate free mode 
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k = aib/v 


K 

I 

m 

n 

<3r 

Or 

T 

V 


moment of inertia of attached mass, relative to elastic axis 

moment of inertia of attached mass, relative to its center of 
gravity 


constant 

number of elastic degrees of freedom; semispan 
attached mass 

number of degrees of freedom 
amplitude of generalized force 
generalized force corresponding to 3{^ 

amplitude of oscillating torque 
free-stream velocity; flutter speed 
flutter speed with zero attached mass 


^r = 


¥ 




X 


y 


->v 


— 

y^ 


El 

GJ 


virtual work 

direct mode xj^j 

conjugate mode ^y^, . • y^] 

direct characteristic mode (nonnormalized) 

conjtigate cheiracteristic mode (nonnormalized) 

displacement of a point of the system- 

displacement of reference point on mass, positive in a 
direction opposite to x 

tending stiffness of -wing 

torsional stiffness of wing 
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p = • • •# ^Iy geometric parameters defining x 

6 torsional deflection of wing, positive nose up 

9 ' angular deflection of body, positive opposite to 9 


X = (cuy/oi)^ 


Xy characteristic number 

V dummy index 

^ spanwlse coordinate of wing 

scalar coefficient of direct characteristic vector . 

Hj. scalar coefficient of conjugate characteristic vector 

oj frequency at flutter 

0 ^ frequency of spring-mounted mass 

bending frequency 
0 J]» reference frequency 

torsional frequency 


HARMONIC RESPONSE FUNCTION 


In accordance with practice, the deflection of a system from some 
initial position of equilibrium is described by n coordinates x^^ 

X 2 , • • Some of these symbols may represent rigid-body displace- 

ments of the system of nonelastic displacements, such as free control- 
Btirface motions. These are called free-body displacements. The others. 
Involving elastic deformation, are called elastic displacements. 

The set of displacements x^, X 2 # • - % arranged in any 

specific way may be conveniently regarded ^ )ctor. The theory of 

reference 1 is essentieilly a geometrical theoi^ based on the vector 
concept of displacements and forces. 
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In accordance with classical mechanics, the selection of coordinates 
defines a system of generalized forces Q 2 ," • • so that the 

work performed in moving these forces through a virtual displacement 
^1 = yi^ Xg = 72^ • • = 7n given by 

W = y^Qi + 7202 + - . . + y^Qn (l) 

The symbol on the right is conveniently regarded as an "inner product" 
of the two vectors y = ^y^, 72 ^ • • and 4 Q 2 , • • 

and written (y, 

Now suppose the system, here an airplane, is subjected to a system 
of harmonic forces and let the symbols q^, q 2 ^ • • In represent 

the amplitudes of the corresponding generalized forces and let Xj, 

X 2 , . . .j Xj^ represent the corresponding an 5 )lttudes of oscillation. 

Then the equations of motion take the form 

y 8j.gXg - oij, ^ bpgXg + q^ = 0 r = 1, 2, . . ., n (2) 

s s 

where the matrix 03^(ars) is the matrix of aero dynamic -plus -inertia 
forces, is a reference frequency, and ~%-^(i’rs) matrix of 

elastic forces. The coefficients a^g depend on the reduced frequency 

k = where b is a reference length (usually the semichord) and v 
the flight speed; the coefficients aj.g are obtained from aerodynamic 
analysis of oscillating flow. 

Dividing equation (2) by cd^ and replacing by X and 

-qp/cu^ by Wj., one has the set of equations 

^ ®Ts^s ~ ^ a^s ~ ^r r = l, 2, ...,n (3) 
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When the external forces are zero, then Wj; = 0 for r = 1, 2, . . n 
and the nonhomogenous set (3) is replaced hy 

H ^s^s - ^ H^rs^s = 0 W 

s s 

The characteristic values Xy associated with equation (4) are the 
zeros of the characteristic determinant ja^g - Xb^gj. To each such zero 
value of Xy there are associated two characteristic vectors or modes, 
one, j called the direct characteristic mode, and 

the other, = (yq'^, • • called the conjugate characteristic 

mode . The direct mode is a solution of equation (4) with X = Xy 

H^s^s'' - \)ZI^rs^s^ = 0 r = 1, 2, . . ., n (5a) 

s s 

whereas the conjugate mode satisfies a similar set of equations 

XIyr''®Ts " '^Y~ yr^'^^rs = 0 8 = 1, 2, . . ., n (5b) 

r r 

The vectors xi'^ and y'^ are not defined uniquely because any 
scalar multiple of these vectors would also satisfy equations (5a) and 
(5b), respectively. In the response function introduced below, a 
"normalization factor" appears which makes it immaterial which solutions 
of equations (5a) and (5b) are used. This factor is 

B(y^ xV) =y~ y^^b^„Xg^ (6) 

r,s 

and represents the work done by the elastic forces set up in the defor- 
mation ^ when these forces are subjected to the displacement 

When aj.g = Sgj. and bj.g = bg^,, the system is conservative and 

the vectors and may be taken as identical. Aerodynamic systems 

are nonconservative in general^ the difference between the vectors 
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and when the two vectors are made coD 5 )arahle hy choosing the same 

component to be unity. Is therefore associated with the energy-absorbing 
or energy-producing characteristic of the system. 

In case is positive, represents the deflection of the 

system when vibrating at the frequency to = When Xy is not 

positive, the mode is sometimes regarded as a similar deflection 

of the system, modified by a suitable dancing coefficient. 

One may, on the basis of equation (5b), regard the conjugate char- 
acteristic mode as the direct characteristic mode of a hypothetical 
system with coefficient matrices aj-g' and l>rs' which are the adjoints 

of the original matrices a^g and bj.g, respectively. Another inter- 
pretation of the conjiigate modes will be made below. 

The harmonic response function is the deflection of the system 
vmder the action of harmonic external forces represented by W]_, Wg, 

. . ., Wjj^. The response function may then be regarded as the deflec- 
tion of the airplane in flight under the action of vibrators. The use 
of the response-function concept has become increasingly In^portant in 
recent years both in flight-testing, as is done for classical dynamic 
stability, and in theory (reference 3) • 

The harmonic response function may be expressed as the stun of 
(a) the response of the system as a free body tuader the action of the 
externeil forces and (b) the sinn of the elastic responses in each degree 

of freedom. The result derived in reference 1 states^ that the harmonic 

, 

response x is 

v=l (xv - 

Here the Integer I represents the number of distinct elastic modes,* 
that is, n - Z is the number of independent free-body degrees of 

freedom. The vector x^'*’^ is the free-body response referred to in 


^For the benefit of the reader who does not have access to refer- 
ence 1, a derivation of equation (7) bas been included in appendix A. 
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Item (a) above; the terms in the summation are the direct characteristic 
vectors^ each appearing with a scalar coefficient 


(Xv - x)B(y'',x^ ) 


r 

(Xv - ^)Y~ yArsV 

r,s 


(8) 


The denominator 
the frequency of the 


of this coefficient contains a term depending on 
forcing oscillation, that is. 


Xv “ X = Xv 


1 

X^Xoi / 


nn^ is of the form of a magnification factor normally appearing in the 
theory of vibration in one degree of freedom. 

The numerator in equation ( 8 ) represents the work fed by the 
external forces into the conjugate mode. Here lies the importance of 

the conjvigate mode; the intensity with which the direct mode appears 

in the response is proportional to the work absorbed by the co njugate 
EodeT" not by the direct mode. This shows that the true nature of the 
couplings due to additional degrees of freedom is to be found from the 
relation of these degrees to the conjugate modes of the original system. 


ADDITION OF A POINT MASS 


Consider an airplane flying at speed v. Suppose that a conplete 
flutter analysis has been made so that the following data are available: 

The characteristic nvonbers X^, the direct modes and the conjugate 

modes 

Now suppose that at a point P of the airplane a concentrated, 

harmonically oscillating force Fe^ is applied. Let x represent 
the displacement of the point P in the direction of the force F. 

Then, within the range of linearity,- 


X = P]_Xj + ^ 2 ^ * 


• + 
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where 3n eviitable geometric constants. The rela- 

tion may he conveniently written as an Inner product of the vector 
^ 32^ • • •} 3n^ and the vector x = X 2 , • • • > 

■ ( 9 ) 

Let ^ qij^ represent the vector with components 

corresponding to the generalized forces developed by the exciting 
force F. In order to obtain the components of consider the work 
developed by the force F under each of the virtual displacements 6x^, 

bxg, . . 5x^. Each such virtual displacement 6Xj. results In a 

displacement Pj,6Xj. of the point P, according to equation (9)j there- 
fore the force F does work F3j.6Xj.. Hence^ 



qr = F3r 


q = FP 



The harmonic response, equation (7), yields for the deflection x 


X 


F_ ^ (y~^^3)(p,x^) _ ^Z+l\ 

0)^ v=l (\i - ^)B(y'',^'') 0)2 ^ ^ 


( 10 ) 


The force F may, for Instance, be the Interaction force produced 
by a mass mounted on the airplane at P. In order to make the results 
more genereil, suppose that the mass m be connected to the airplane at 
P by a spring (fig. l(a)). Let 0)^3 represent the frequency of the 

mass on the spring, P held fixed. Then the deflection x' of the 
end of the spring under the reaction F Is 



X 


m 


( 11 ) 
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Since X = -x', con 5 >arlson of equations (lO) and (U) gives 


lL-1- + jJ\ = JL 

03^ 0^2 j ^2 ^ _ x^b('^'',xV) ^ / 


m 


Dividing throu^ by F/(j^, 


i/ 1 + ^ \ _ (y^ }^)(^^) JLfig' 

A 0^^ I ^ hi ^ 


(12a) 


In the case where there are no free modes^ 2 = n, and equation (l2a) 
is replaced by 




^ (Xv - X)B(^V,-iV) 


(12b) 


In case the mass is rigidly attached to the point P, is 

infinite and the term is zero. 

Consider the case in which the mass is attached with a spring to 
the wing. Suppose, in addition, that the deflection of the wing is 
described by two degrees of freedom, one bending and represented by 
fjj^(l), and one torsion represented by fo,(|). Then if e^ represents 

the chordwise coordinate of the point P relative to the elastic axis, 
positive rearward 


t = (fjid), (13) 

The right sides of equations (12) are functions only of v/oib and 
(£)/(Dj.. Hence equation (12) may be solved, in principle, for m and 

for each vaJLue of v/oto. In practice, it appears desirable to 

carry out the solution by a graphical construction. In this procedure, 
one regards the right side of equation (12), for each value of v/oib. 
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as a function of tlie real parameter X = The function Is plotted 

In the complex plane and the abscissa of Intersection of the resulting 
plot on the axis of reals Is equal to 


iL . . i(.x . i 

m\ m\ X 


(14) 


for -which flu-tter -will occur. The parameter X at the Intersection 
defines the frequency of flutter and the -value of m corresponding to 
equation (l4) the mass required to maintain flutter. As In conventional 
practice (reference 4), the mass Is to be Interpreted as the nondlmen- 
slonal mass m/jtpb3. 

The graphical construction la facilitated by noting that each of 
-the summands In equation (12) has the form c/(Xi - X) where c Is a 
constant and X^ Is a characteristic value. Both c and X^ = a]_ + Ibj 
are conq)lex numbers In general. A plot of the fraction c/(Xi - X), with 
X a real parameter -varying from min us Inflnl-ty to plus Inflnl-by, yields 
a circle passing through the origin of the coiiq)lex plane. To construct 
the circle In the complex plane: 

(a) Form the diametral -vector c/lb and draw the corresponding 
circle through -the origin (fig. 2). The diametral vector represents 
the complex number of largest modultis given by c/(ai + Ibx - X), 
attained when X = a^. 

(b) Give a sense to -fche direction of motion around the circle as 
X -varies from minus Infinity to plus Infinity according to the rule: 

If bx > 0, the circle Is traversed clockwise 

If bx < Oj the circle Is traversed counterclockwise 

(c) Calibrate the circle. Eeferrlng to figure 2^ use any convenient 
scale^ lay out distance |bx| along the diametral -vector, and draw 
line PP perpendicular to -the diametral -vector. Give the same sense to 
points on PP as has been established on the circle. To point Q, attach 
the -value X = a^. To any other point Q' on line PP, a-ttach the 

-value X = a^ ± iQQ'l, where jOQ'l Is the distance from Q to Q' 
measured to the scale chosen abo-ve for |T3 x|^ and the plus or minus 
sign Is selected to correspond -with the sense of increasing or decreasing 
-values of X. Draw the radial -vector from the origin through Q’ . The 
•value of X at Q’ is then attached to the intersection R of the radlEil 
vector on the circle. 



NACA TN 251^0 


13 


The procedure is carried out for each constituent of the response 
function. The resulting vectors OR, corresponding to the same values 
of X, are added vectorial 1y.^ The locus of the sum in the complex 
plane is a plot of the harmonic response which appears bn the right side 
of equation (l2b). 

The analysis described above is illustrated in appendix B. 


ABDITIOR OF A DISTRIBUTED MASS 


In the preceding section, the mass was regarded as concentrated in 
a point; that is, the moment of inertia I of the mass relative to its 
center of gravity is zero. In this section, ttie equation corresponding 
to equation (l2) but for a distributed mass (l ^ O), rigidly attached to 
the original system, will be derived. 

There are now two conditions required that the deflection of the 
mass and the system be compatible. Two equations of compatibility are 
therefore derived . 

If the original system is described only by two degrees of freedom, 
then the analysis of the mass effect as described below is not a con- 
venient device because the conditions of compatibility are equal in 
number to the original equations of flutter. On the other hand, if a 
large number of degrees of freedom were required to describe the original 
system, then the numerical work required to treat the two equations of 
compatibility may be significantly less than a reworking of the entire 
flutter analysis to include the mass effect. 

In order to make. the analysis definite, suppose the mass is attached 
to an airplane wing which can vibrate only transverse to the plane of 
the wing. • Let P be a point common to the body and mass (see fig. 1(b)) 
and 

X vertical displacement of wing at P 

x' vertical displacement of body at P 

6 torsional displacement of wing at P 

0' torsional displacement of body at P 

The sign convention is such that x' and 0’ are positive in directions 
opposite to X and 0. 

2ln order to perform the vectorlcal addition graphically, it is 
desirable that the individual circles, and therefore the diametral vectors, 
be drawn to the same scale. 
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There are nov geometric parameters an^ * 3 ^ such that 3 


= (1^4) 
= ( 1 ^, 4 ) 


At point Pj apply a harmonic force F and torque T. Let 

vertical displacements of wing at point P due to F and T, 
respectively 

torsional displacements of wing due to F and T, respectively 

vertical displacements of body at point P due to reversed F 
and T, respectively 

't torsional displacements of body due to reversed F and T, 
respectively 

Under the action of F, the forcing vector v is w' = — — P . 

Under the action of T, the forcing vectpr is Therefore, 

or 


^ = .JLy- *1 ^Z+l) 


^ ■ 'o)2 ^ (Xv - x)B(rA^) ^ 


= -X. _ F 1^2 ^z+i') 

0)2 (\y - ^ 


3 ln the particular case considered here, of wing bending and wing 
torsion. 
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e 


T 


jL 'v X(^2 -^z+i) 


_1_ _I_ 
03 ^ ml 


F 


e-p' 



(16) 


x^’ 




The condition of compatihility results in. a pair of simultaneous 
equations 


Xp + Xp’ + X^ + Xqi' = 0 

0p + 0p ' + 0^ ' + drjt — 0 


(IT) 


These equations can he expressed in terms of the characteristics 
of the wing and mass by substituting from equations (15) and (16) . 
Restricting the discussion, for slii^jliclty, to the case where there are 
no ffee modes present (2 = n), one has 


0 = F 


0 = F 



>(18) 


1 

I 


where the factor l/cu^ has been canceled out. 
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Let A represent the determinant of the system (I8): 


Ali 

^ ®ii 


Ai 2 

Bi 2 

~r 

•• "Xf •* ^ 

_ 

ml 

“T 

” X X^ “ X 

A21 

B21 


A22 

®22 

A«j|^ “* Xi ” "X* 

+ — 
I 

T 

Xj^ “ X X 


where 


Then the condition for the existence of a nontrivial solution of equa- 
tion (18) reads 


A = 0 


(19) 


It will now he assumed that the matrix b^,g is the unit matrix 
and that the direct and conjugate modes have been normalized so that 


- 1 


V = 1,2 


( 20 ) 


These assumptions represent no additional limitation. 
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In view of equation (l 8 )j the expansion of the determinant A 
becomes 


A = 


All 

A21 


Ai2 

A22 


(Xl - X)2 + 


®11 

®21 


®12 

®22 


j {\2 - X )2 


+ 



B 


12 




The first three terms may be simplified as follows: 


All A]^ 

A21 A22 






0 

H 



0 ( 1 ^,?^) 

X 



0 


in view of the vanishing of the second determinant in the product. 
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Similarly |Bi.s| = 0. Also, 

All ®12 

(|i,^) {t\P) (p-,P) 0 

^ (P,^) 0 (^,p2) 

0 

(T,^) (|2,#) 0 

(p-,^ R^,!^ 0 (f,P) 0 

p,-^) p,1^) |_ 0 (^,1^ ■ 0 (|l,p2)_ 

(?-,?-) ^ (P,P) 

p,n {^^^) "" (p,p) i^,-^^) 

^x Pa^ ^ ^ ^x ^2^ 

XXX 
^X ^ ^ ^2^ ^2^ 

From the orthogonality property (appendix A) and the normalization 
(equation (20)) and in view of the fact that matrix B has "been assumed 
to he the unit matrix, 

b(P-,'x^) = = 1 

b(^,^) = = 1 

b(^,^) = (■^,x2) = 0 

b(^,^) = = 0 
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Therefore, 


Hence, 


x'l X 1 


7 l^ 


(&,F) 


X 2 y 2 
Xi x^ 

X 

Y 2 ^ Y 2 ^ 





1 0 
0 1 


= 1 


■^1 ®12 


®u -^12 




+ 


= 


^21 ®22 


®21 -^2 


Pi" P2" 


( 24 ) 


Finally, 


._L 

ml 

be 

I 


be 

I 

_ 1 

f 


-2 " -2 
ml I 




I - 0,e)2] = J_ 

- . J. Im 


Combining equations ( 2 l) to ( 24 ) and multiplying through by I, the 
equation of flutter (19) becomes , • 


lA = I 



^11 


(Xg - x) + i = 0 


be 




/(^l "^) 


(25) 
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The first term on the right may he separated into partial fractions 
and combined "with the second and third. Equation (25)^ so modified, has 
the same form as equation (l2b) for the case in which is Infinite 

(rigid attachment) . Therefore the graphical solution of equation (l2b) 
is directly applicable to equation ( 25 ). 

Such a graphical solution is illustrated, for the numerical example 
discussed above, in figiire 2 (h). 


RUMEKECAL ILLIETRATION 


The method outlined has heen applied to obtain the flutter speeds 
under figure 4 of reference 2. Here the weight is attached at 
= 0.500 for various span positions. A detailed calculation is 
presented in appendix B. The results of an analysis are presented in 
figures 3 to 7* 

The specific mass addition with which this report is concerned, 

that is, — ? = 357^ is represented by a horizontal line in figure 3 . 

jrpb^ 

The intersection at each span position gives the corresponding value 

of ~ at which flutter will occur. For these values of — the 
oh a)b 

values of \ are deteormined from figure 4. Table 1 contains the values 

Y 

of ^ and X read from these curves. 

Figure 5 "was plotted in a form con^jarable with figure 4 of refer- 
ence 2. However, the definition of Vq differs from that of refer- 
ence 1. Here Vq is the flutter speed of the wing without attached 
mass according to theoiy or test as the case may be. The cuorves which 

are being conq)ared are (-^) and {— ) 

' 0 / theory '^ 0 / test 

Figure 6 presents a comparison of the reduced flutter speeds 
obtained by the method of this report and reference 2. Figure 7 shows 
a comparison of frequency ratios. In particular it is shown that the 
effect of I, in this case, is small. 

The effect of flexibility of the suspension is considered in 
appendix C. 


Purdue University 

Leifayette, Ind., July 20, 1950 
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APPENDIX A 

DERIVATION OF HARMONIC RESPONSE FUNCTION 


In the following paragraphs, a derivation of equation (7) is 
presented. This derivation is based upon a number of assun 5 )tions which 
are discussed in reference 1. They are: 

(a) Let e and f represent direct and conjugate characteristic 
vectors associated with a given characteristic value. Then,^ 

B(f,t) / 0 

(b) Let g and h be free-body direct and conjugate vector s^ such 

that 

^ ^rs^s “ ® 
s 

hpbj.a = 0 

r 

— > 

Then, for each vector g, there is at least one vector h such that 

H l^s^sgs ^ 0 

r,s 

(c) The characteristic determinant |aj.g - Xbj,g| does not vanish 
Identically. 

Consider the characteristic-value problem 

^ ®rs^s “ ^ ^ ^ra^R ~ ® r = 1, 2, . . ., n (Al) 

8 S 


^*This condition teikes the place of "positive-definiteness" of the 
energy, required in the emalysis of conservative systems. The normal- 
ization condition associated with the vectors "e and T is introduced 
below. 

Sectors defining motions with zero elastic potential energy. 
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Equations (Al) represent the balancing between the Inertla-plus- 
aerodynamlc forces represented by the matrix a^g and the elastic 
forces, represented by -X(brs) developed In the mode 
X = xj, X2, . . x^. 

Another Interpretation of equations (Al) Is obtained by multi- 
plying the rth equation by y^ and summing on r_. The vector 
y = yij 72.3 • * ' 3 - y XL represents a virtual displacement. Introducing 
the notation. 


A(y,x) =y~ Yr^-s^s 
rs 

B(y;,x) = yr^rs^s 
rs 




(A2) 




one may write the result In the form 


A(y,x) - XB(y,x) = 0 (A3) 

for al l values of y. In this form, the equation of dynamic equlllb- 
rliim states that the sxmi of the works performed by the Inertla-plus- 
aerodynamlc forces and elastic forces, when displaced- over an axlbltrary 
virtual displacement, is zero (d’Alembert’s principle).. 

In the special case where bj.g Is the unit matrix. 


B(y,x) = ^ y^x^ 
r=l 


Then B(y,x) will be written simply as (^,^ . In this case, 

B(yjx) = (y,x) has the Interpretation as an "inner product” between 
two vectors y and In two and three dimensions, (y^x) Is then 

equal to the product of the lengths of 5c and y by the cosine of 

the Included angle. In the case of arbitrary values of bj-s^ ^ similar 

Interpretation will be made; that is, B(y,x) will be regarded as an 
Inner .product between the conjugate mode and the direct mode x. 
Further, “y and "x will be orthogonal with- regard to b^s if 


(Alf) 


B(J,2) . 0 
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It is to be emphasized thatj in using this concept of an inner product^ 
the direct vector must appear to the second argument nni^ the conjugate 
vector in the first argument. 

Similarly _^(y,x) defines an Inner product. It, is conceivable 
that a vector y may be orthogonal to with respect to bj.g but ' 

not orthogonal with respect to aj.g and vice versa. But equation (A3) 

states that^ if x = is a direct characteristic mode corresponding 
to a characteristic root X = ' then 

A(yje^) = X^B(y,el) (A5) 

for all values of “y. Therefore, if is orthogonal to "el with 
respect to b^g, then the right side of equation (A5) is zero and 

so also is the left) then y is orthogonal to also with respect 

to a^g. 

One may, along with equations (Al), consider the set 


y ^r^rs “ ^ y^r^rs ~ ® s - 1, 2, . . ., n (a6) 

r r 

Multiplying the sth eqmtion by Xg . and summipg on s, one has again 
equation (A3) which must be satisfied for arbitrary, values of. x. 

Since X = is a characteristic root of equation (a 6), it follows 

that there is a nontrivial solution y = f 1 of equation (a 6) for 
X = Xj^. Therefore, analogous to equation (A5), 

a(?-,x) = Xj^b(P-,1^ (A7) 


for all values of "x. 

Whereas to the classiceQ. theory of conservative systems the 

vectors and may be taken as identical, the nonconservative 

character of aerodynamic forces requires that ~e^ and ¥*■ be treated 
separately. 

Let ’jT be an arbitrary vector representing a direct mode. Then 
resolve x into components "parallel" and orthogonal with reference 
to b^g, to That is, decompose IT to the form 

S = + S<2) (AS) 


T 
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■wiiere | is a scalar and is sucli that 

= 0 (A9) 

Such a deccanposition. is possible because equation (A9) is equivalent to 


- 1 ,:^) = 0 

B(r*-,x) - = 0 

= b(?i,1?-)/b(?i,?i) 


(AlO) 


At the same time, one may resolve y into coniponents, parallel 


and orthogonal, "wi"th reference to b^g, to e^. 


y = Tijfl 

(All) 

= 0 

(A12) 

Til = B(?,?l)/B(f3-,tl) 

(A13) 


Since ^ 0 (Assumption (a)), one may assume, for con- 

venience, that 

B(?-,e^) = 1 (A14) 


a condition ("normalization condition") "which may be met by multiplica- 
tion of "e^ or f^ by a suitably chosen scalar. 


The orthogonali"ty relations (A9) and (A12) "with respect to b^g 
inqoly orthogonality "with respect to apg, as stated abo"ve. 


A(fl,^2)) = 0 


(A15) 
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(A16) 


so ’that there is zero coupling between the first characteristic inodes 
and the orthogonal modes. 

In view of equations (A9), (A12), (A15), and (AI 6 ), and ?■ 

are particularly well-suited for this construction of a set of coordinate 
vectors^ one for the direct modes and one for the conjugate modes, with 
reference to which A(y,x) and B(y,5) take simple forms. For, 

= ni^lA(?'-,el) + A(y(2),l(2)j 
Similarly, using equation (Al4) 

B(J,S) = + E(y(2),t(2)) 

= I 1 S 1 + B^(2),x(2)) (A18) 

Substituting y = ?*- in equation (A5) and using equation (A14), 

a(|i,^) = (A 19 ) 


and equation (Al?) takes the form 

A(y,x) = XiTi2^52. + A(y(2)^x(2)) (A20) 

In virtue of the orthogonality conditions (A9) and-(A12), the 

■*( 2 ) -^ 2 ) 

vectors x' and each have n - 1 degrees of freedom. Every 

value of x^“^^ may then be represented as a linear combination of 



n - 1 vectors Similarly, every value of y(^) is a linear com- ^ 

bination of vectors y^, • • Then an arbitrary direct mode x may be 

vritten 


X =■ + . . . + Ir. „ 

JU CJ- c^u-x 


and an arbitrary mode y 


j = ih 


51 , _ ^=21 


t • • 


•♦P.T1-1 

’I2,n-1^“'“ “ 


With the aid of equations (A18) and (A20), equation (A3) may now be written 

^l^lSl + ^ s 

r,B-l ' [_ r, 8=1 * _ 


= 0 (A21) 


nitting tne coefficients of tj^ and t |2 j. equal to zero, a set of linear equations 
in and Ig^j. is obtained. The characteristic equation is now 


X.3^ - 0 0 

0 -xb^A,-^A) ^ 

... 0 
. . . 

0 ... 

. . . 


1 
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It can nov "be seen that the characteristic roots of the original 
problem in n degrees are identical with together with the set 

of the (n. - l) -degree problem corresponding to 

Since the asstnnptions (a) to (c), assumed valid for the n-degree 
problem, hold for the (n - l) -degree problem, one may continue the 
reduction. Uamely, to a characteristic root X^ of the characteristic 

determinant associated with there is . 

a direct characteristic mode and a conjugate characteristic mode 

Assuming that they are normalized. 


= 1 


then it is possible to resolve 

B^2^'5(3)j = 0 
{2 = b(?2^-5(2)) 

At the same time, 

b(5(2>3(2)) 


x(^) and y^^^ in the forme 

B^y(3)^-g2^ = 0 

^2 = 

= ^2^2^ A^y(3),?(3)) 

= ^2I2+b(^3),^(3)) 


The characteristic equation then takes the form 


0 = 


X*j^ — \ 
0 


0 


0 


0 I 0 
X^2 “X# I 0 
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One may continue tlie reduction process until a.1 1 the roots of the 
characteristic equation have been used. If the reduction has been 
made n times, then. 



A(y,t) = Mrir 

r=l 

. tl \lr 

1 

and the characteristic equation is 


1 

0 

0 

X 2 “ ^ 

... 0 

... 0 

0 

0 

• • • ^ 


On the other hand, if eifter Z reductions the characteristic 
determinant is a constant 0, by assumption (c)), then it can be 
shown (reference l) that the reduced bpg is the zero matrix. Hence, 


X 


21 

1 




s 

A(^,^ = ^ Mrir A^y^+l,x^+l) 

B(y,x) =YlTir5r 

r=l 


> (A22) 



NACA TN 2540 


29 


The orthogonality 


conditions 


6 


prevail; 


A(fr,ts) = 

= 0 



= b(P,$^+i) . 0 

r = 


= = 0 

s = 


The characteristic equation is. 


2 , 

2 , 


I 


K (A23) 


I 


j 



where the terms indicated with asterisks are independent of 

From the last set of equations (A23), it follows that 
and y^"^^ are free-hody modes. 

With the reduction shown in equation (A22), it becomes a sln 5 )le 
matter to solve equations (3) which are more conveniently treated in 
the form of an equation of virtual work 

AtJ,?) - (A24) 

for all values of y. 

°The orthogonality described here is sometimes referred to as "bi- 
orthogonallty. " For Instance, see reference 5 where there is presented 
a method of analysis, similar to that in reference 1, for calculating 
the effects on flutter of small changes in the parameters. 
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Siibstl-tirtlrig from equation (A22), 

*xVlr + ^ Irlr = (ZI + y*+l,v) 

for all •values of and an values of Hence, 

ih: - !^)Sr = (A25) 

and, for all -values of 

a(P+1^5Z+1) . (J2+1,?) (A 26 ) 

From equation (A25)j 

= (^,v)/(X2 - ^ (A27) 

Equation (A26) defines the free-hody motion (elastic degrees regarded 
as frozen) \mder the action of the external forces and may be reduced 
to a set of linear equations. Substituting from equation (A27) to 
equation (A22), 

1 = y ty + (a 28 ) 

Z ^ - X. 

In case the normalization condition (A14) does not hold, then it 
can be shown that eqmtion (A28) is modified as shown in equation (7)» 
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APPENDIX B 
SAMPLE CALCULATION 


In "this appendix^ a sample calculation Is canried out for the 
example discussed in reference 2. G?he case considered corresponds to 
e^ = 0.500 and cUq = « (rigid attachment). 


Homogenous Solution 

The flutter determinant to be solved in two degrees of freedom, 
according to the notation of reference 3, pages 62 to 63, is 


A B 


an - Xbn 

®-12 

D E 


a^ 

0^22 "" 


= 0 


For a wing with the 


following properties: 


Chord, ft . . . 
Length, ft . . 
Aspect ratio . 
Taper ratio . 
Airfoil section 
7, slugs/ft 
I^ slug-ft2/ft 

El, lb-ft2 . . , 
GJ, lb-ft2 . . , 
m, slugs . . . , 
I, slugs . . . , 
radians/ sec 
ci^ radians/ sec 



1 

NACA 16-010 
. 0.02702 
. 0.00080 
-0.126 
977.08 
480.56 
0.0988 
. 0.00452 
4l.l 
299.6 


The shapes of the deflection curves in bending f^(|) and torsion 
fo^(^) along the span are 

= i^cosh 1.87 Y - cos 1.87 ij+ 0.734^sin I.87 i - slnh I.87 

fjl) = Bln II 



so that 


U) 

ro 









fn(6)fa^6)'^& 


1.4-114 


.. p 

V n 

Consider as an exainple ^ = 8.33. Define >. = — The values of the terms in the 
ae*fce2Tnlriairfc are 


^11 " ^^11 

^12 


( 91.6464 - 40 . 31551 ) - X(l. 841762) 

(-476.7902 + 34.89841) 

&-JL 

Q.rv^ 

tzc . ££ 


(9.97095 + 21=2750541) 

^ klT 51 R _ to R*^^or«^ ^ \ f AllA)lll' 7 ^ 

\ .jww* 1* 1,^.^ - / “ r*\ j 


Solving a quadratic equation for the tasic characteristic value X, the values obtained for the 
two modes are . 

\ = 5.9082 + 1.0438i 


Xg .. 49.6991 - 24.2515i 


NACA TN 2:540 
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Enoiigh information is now available to con^pute the direct and 
conjugate characteristic modes from the vlrtrual -work equations thus: 


6 Work = 0 = 


0 = 

&21 




22 


^1 

Xb22)\^. 


Let 

y^V = x^V = 1 

Then, 


X V = _ (^11 - Vn) 

2 . , ^12 



( a^i - Xbii ) 

a2i 


^21 

( a22 - ^^22) 

^12 

(a22 “ ^^22) 


Table 2 has been- conq>iled of all the characteristic modes thus 
conq)uted . 


Besponse Fxmction - Mass Addition 


The response function developed in appendix A when e3q)anded is 

B(yv,xv). 

m/npb^ 

[yi^hU) + y2^Vg(S)] [xi^^h(^) + ^ 

(Xl - X)B(yl,xl) , 



[y2^h(^) + y2^ewfg(^)l + ^^^yfc^Cg j] 

(^2 - x)B(y^,x2) 



Since y-j^v m x^V = 

R.F „ [^h(0 + [fii(^) + [^(6) + ^^eyfjg)] 

(Xl - ^)B(yl,xl) (Xg - X)B(y2,x2) 

% ^ Kg 

^ ^*2 "" ^ 

Consider again ^ = 8*33 and a span position of 75 percent (| ■=■ 3')j then Kj_ anfl K2 are 
coB^juted from the following data: 

fj^(3) » 0.6550^^5 

e^j3) ° 0.461940 

Substituting, 


[0.65504^ + (0.172 + 3»868i)(o.46i94og [0.695045 + (0.175 - O.Q7^7i) (0.46194 j 

(Xl - x)(i 8.6357 + 34.48831) 


To. 655045 + (-0.171 - 0.07401) (0. 461940 )1 To. 655045 + (-0.000405 + 0.009301) (0.46l94o)~ 

■ K ■ ■ ■ ■ .1 I. .11 . . , , , ■ ^ h-T _ _ -J 

(Xg - x)( 1-88 i 7 - o.o8lo4l) 

0.0362119 + 0.002116541 0.200690 - 0.001938561 

Xi - X X2 - X 


K -^ = 0.0362119 + 0.002116541 
Kg = 0.200690 - 0.001938561 
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The values of the constants and are shown in table 3* 

Having found K-i and the values of X and ' — p — =■ are 

^ m/jrpb^ 

determined as shown graphically in figure 2. The data calculated in 
this manner are conqjiled in table 1. For these con 5 >utations natural 
frequencies and found in test and corrected for the effect 

of apparent mass have been used (table h) . 
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APPENDIX C 

FLEXIBILITY OF MASS SUSPENSION 


The numerical example^ figures 1 and 3 to 6, has been carried out 
for the case in which the mass m is rigidly attached to the wing. 

The results plotted in figures 3 and h are still valid if one replaces 
m hy 


m 




1 , 

1 

X 



Then, for a given value of m, the data of figures 1 and 3 niay be cross- 
plotted to give a relation between the flutter speed v. 

This has been done for the case where the mass is attached at 
50-percent span and the results are presented in figure 8. When the 
stiffness of the mass suspension decreases, the flutter speed decreases. 
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TABLE 2 


DUfflCT ABD CONJUGATE MODES 


v/a* 


X V 


XgV 

yi^ 

ygv 

B(yv,xv) 

5-00 

\ 

2.8155 - 0.14841 

1,0 

0.583 - 0.1421 

1.0 

- 2.440 + 6.0551 

-27.5064 + 201.86021 

Xq 

51.3996 - 12.27451 

1.0 

0.00203 + 0.005031 

1.0 

-0.0573 - 0.01401 

1.83937 - 0.0164921 

6.25 

h 

3.8861 + 0.o8q94i 

1.0 

0.345 - 0.09941 

1.0 

-1.099 + 5.0761 

8.3356 + 96.90191 

Xg 

50.5943 - 16.43811 

1.0 



0 . 00144 - + 0.006661 

1.0 

- 0.0946 - 0.02731 

1.8441 - 0.034891 

8.33 


5,9882 + 1.04381 

1.0 

0.175 - 0 . 0757 i 

1.0 

0.172 -f 3.8681 

18,6357 + 34.48831 

Xg 

49.6991 - 24.25151 

1.0 

-o.ooo 4 o 5 + 0.009301 

1 . 0 , 

-0.171 - 0.07401 

1.8812 - 0.08l04l 

10.00 


717401 + 2.35621 

1.0 

0.113 - 0.06711 

1.0 

0.764 + 3.1961 

17.4791 + 16,13151 

Xg 

49.5819 - 31.23091 

1.0 

-0.00250 + 0 . 01051 ' 

1.0 

-0.231 - 0.1381 

1.9470 - 0.10821 

12.50 


9.7604 + 4.73691 

1.0 

0.0618 - 0.05801 

1.0 

1.296 + 2.5i4i 

\ 

13.5918 + 4.15641 

Xg 

51.0777 “ 42.25721 

1.0 

-0,00573 + 0.01111 

1.0 

-0.302 - 0.2621 

2.0834 - 0.096711 

16.67 

^1 

11.6265 + 8.77351 

1.0 

0.035 7 0.04691 

1.0 

1.737 + 1.8501 

9.5504 - 8.78151 

^ . 

56.7962 - 61.02201 

1.0 

-0.00954 + 0.01021 

1.0 

-0.361 - 0.4841 

2.2775 + 0.049281 




LO 

VD 
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1 


25-percent 




Kg 

5.00 

0.00115529. + 
0.0000685931 

0.00456109 - 
0.0000562521 

6.25 

0.00164108 + 
0 . 0001412841 

o.oo 4 l 48 l 6 - 
0.0001420871 

0.33 

' 

0.00245659 + 

r\ r\ry^r^r>/C/^Fr!i-f 1 

V * 1 

0.00332996 - 
0 • 0005238651 

— ■■ 1 

10.00 1 

0.00310723 + 1 
0.001096581- 

0.00268993 - 

0.0010241 

12.50 

0.00387492 + 
0 . 0023775*^1 1 

0.00194353 - 

0.002172611 

16.67 ' 

0.000906428 
0.004348551 i 

0.00115146 - 

0.003062651 


table 3 

VALUES OF COKSTAKDS AHD Kg 


% ^ 
X]_ - X ■*■ X2 - X 


50-percent span 

75-percent opan 

100-percent epon 


Kg 

. % 

«2 

' % 

Kg 

0.00591099 + 

0.0588210 - 

0.0128929 - 

0.225415 - 

0.0191386 - 

0.531505 + 

0.0000763931 

0.0000801211 

0.000504001 

0.00109631 

0.00226521 

0.00227061 

0.00850877 + 

0.0559070 - 

0.0196935 - 

0.217386 - 

0.0308311 - 

0.516939 - 

0.000321311 

0.000324951 

0.000655261 

0.0006486921 

0.004100581 

0. 004093451 

0.0144091 + 

Q.o 1<99653 - 

0.0362119 + 

0.200690 - 

0,. 0614892 - 

0.485967 + 

n nrvoi 

ri ruv 3 nf^iv 7 ri 4 

W * \y>./^V./WV/ 1 Wak 

0 . 002116541 

0.001938651 

r» nn -31 kn l 

Ci r^no*^‘ 5 inCH 

0.0196925 + 

0.0447656 - 

0.0523205 + 

. - - 

0.184858 - 

0.0902713 + 

0.454338 - 


0.005470381 

0,009662761. 

0.009770271 

0.01062461 

0.007760971 

0.0267443 + 

0,0380415 - 

0.0753990 + 

0.163083 - 

0. 142436 + 

0.408613 - 

O.Ol 4 l 4991 

0.01273461 

0.03287941 

0.02905991 

0.04826591 

0.04146871 

0.00838717 + 

0.0306178 - 

0.0282219 + 

0.130531 - 

0.0621975 + 

0.356243 - 

0.03011491 

0.02546841 

0,08287861 

0.06585941 

0.1514681 

0.1116231 




i 

J 

f 
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03ABIE 4 

COMPAEISON OF CALCULATED AKD TEST FREQUENCIES 


\ 

Calculated 

Test 

Test (corrected) 
(1) 

First bending (radiajis/sec) 

41.6 

40.5 

41.1 

Second tending (radians/sec) 

262 

246 

249.8 

First torsion (radians/aec) 

304.4 

297.2 

299.6 

Second torsion (radians/sec) 

914 


, 914 


^eat data are corrected for the effect of apparent mass. 







F 

;ri‘buted mass, rigidly attached. 

for relative deflection of mass and wing. 




Vaj Point massj I = o. 

Figure 2.- Kesponse function. ~ - 8 r. 

aib ~ °*33> 75-percent span 

Fractions represent scale of calibration line. 

^ 0:003H6-5ti , 0^00690 - o.ooi9«^f;i 
5.9«S2 + 1.04381 . 1 5979991 - 24.25151 . x ’ — 

my 

^ = 6.0, diam. = 0.00203 - 0.0347i- X - to 7 

dlan. = 0.OCOC8 . 0.0082T51. “fStVo^^il'for clrcl. 






("b) Distributed massj I = 0.00452. 


0.052466 + O.OOT9^51i 0.191216 - O.OOTT9^1i 
5.9882 + l.o438i - \ * 49.6991 - 24.25151 - 1 




/ Jtpb^ 


X = 6.0, diam. = 0.007621 - 0.05028ij X = 49.7, 

diam. = 0.0003214 + O.OO 78851 . 1 Incb = O.OI 8 for circles 


Figure 2.- Concluded. 





Re^c/uce,d ap<ze,cf^ 


(■fa) Distributed massj I = 0.00452. 
Figure 3*” Concluded. 




) 


] 





Rec/uced -fluHer spaec/-, 




ro 


§ 


1 


(a) Point mass; 1=0. 

Figure 4 .- Frequency ratio against reduced speed. 


•f=- 





Reduced epee cub 

(b) Distributed maasj I ■= 0.00452. 

Figure 4.- Concluded. 





F hi ter speed raiio^ 






A. Re^ sponse method CX’‘0) 

S. FesponsQ JT}ethod^O.OC452) 
q. N/\ CA m/90Z, test 
n /VACA TNIBOZ, colcuhfed 


Posit ion^ percerri span 


IT -1 »-«-i 1 + 

X' IXX ^ «/*”" X'JUIXW 


OTt 


I 


o^C^ ill ^^zm 



Reduced fluffs r sfoeed 


Ul 

o 
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Figure 6.- Reduced flutter speed- 


on^s Hi v:jvm 






